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RESUMEN 

Se presentan diagramas de diagnostico de densidad electronica - excitacion 
para una muestra de 613 nebulosas planetarias. En base a esta extensa muestra se 
definen nuevos limites estadi'sticos para la distribucion de nebulosas planetarias en 
Ids pianos log [Ha/[Sii]] ws log [Ha/[Nii]], log [Ha/[Sii]] vs [Sii] AA6717/6731 y log 
[Ha/[Nii]] vs [Sii] AA6717/6731. Los diagramas proveen una buena representacion 
de los intervalos de condiciones fi'sicas, indicados por estos cocientes de Imeas de 
emision, presentes en nebulosas planetarias en diferentes estados de evolucion. 

ABSTRACT 

Diagnostic diagrams of electron density - excitation for a sample of 613 plan- 
etary nebulae are presented. The present extensive sample allows the definition of 
new statistical limits for the distribution of planetary nebulae in the log [Ha/[Sii]] 
vs log [Ha/[Nii]], log [Ha/[Sii]] vs [Sii] AA6717/6731 and log [Ha/[Nil]] vs [Sli] 
AA6717/6731 planes. The diagrams provide a good representation of the ranges 
of physical conditions, indicated by these emission line ratios, present in planetary 
nebulae during different evolutionary stages. 
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1. INTRODUCTION 

Plots of emission line ratios provide useful diag- 
nostics of the physical conditions and identification 
of classes of ionized gaseous nebulae. Electron den- 
sity - excitation diagrams were introduced by Sab- 
badin, Minello and Bianchini (1977), SMB77 here- 
after, with the specific purpose of clarifying the na- 
ture of the nebula S 176. These diagrams compare 
the relative line intensities of the Ha/ [Nil], [Sii] 
AA6717/6731 and Ha/[Sii] ratios observed in super- 
nova remnants, planetary nebulae (PNe) and HII re- 
gions. These nebulae are expected to occupy distinct 
regions in these diagrams as a consequence of the 
different physical processes operating within them. 
For example, supernova remnants are mainly shock 
excited nebulae where low Ha/ [Nil] and Ha/[Sii] 
ratios are expected across a range of electron den- 
sities. H II regions are photoionized nebulae that 
show a restricted range of excitation conditions and 
electron densities in their generally diluted environ- 
ments. Meanwhile, planetary nebulae show a wide 
range of excitation conditions and electron densities 
depending upon their evolutionary stage, the possi- 
ble presence of coUimated outflows and shocks and 



the progenitor mass that may influence the observed 
Ha/ [Nil] ratio. 

Recently, Phillips (2004) has discussed the dis- 
tribution and significance of other relevant emission 
line ratios in PNe. Emission line ratio plots have 
also been used to distinguish planetary from sym- 
biotic nebulae (e.g. Gutierrez Moreno, Moreno & 
Cortez, 1995). The SMB77 diagrams have proven to 
be a useful diagnostic tool for PNe and have been 
used widely in a number of works (e.g. Barral et 
al. 1982; Lopez and Meaburn 1983; Garcia Lario 
et al. 1991; Tajitsu et al. 1999). In the case of 
PNe the data in these diagrams often fall outside 
the original boundaries fit by SMB77, most likley as 
a consquence of the limited sample available then to 
set the diagram limits. The information used orig- 
inally by SMB77 to define the parameter space for 
the PNe was obtained from the limited compilation 
made by Kaler (1976), where relative emission line 
intensities observed in planetary and diffuse nebu- 
lae are listed from different sources. The estimated 
available number of PNe in that catalogue with the 
required information is of only about 40 objects. The 
information for H II regions in SMB77 was also ob- 
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tained from that same source, whereas the data for 
Supernova remnants from Daltabuit et al. (1976), 
Sabbadin and D'Odorico (1976) and Dopita (1976). 

Here we have taken advantage of line intensities 
listed in the Strasbourg Catalogue of Galactic Plan- 
etary Nebulae, part II, SCGPN II hereafter, (Acker 
ct al. 1992) to explore the empirical location of PNe 
in the SMB77 diagrams, but now based upon a much 
larger and homogeneous sample. 

2. THE SAMPLE 

In order to update the parameter intervals of the 
log [Ha/[Nii]], [Sii] AA6717/6731 and log [Ha/[Sii]] 
line ratios for PNe, we have compiled a database with 
613 objects from the SCGPN II that provide the line 
intensities required in the diagrams. 

The Ha/ [Nil] ratio involves both [N II] hues, 
namely A6548 & A6584, whereas the SCGPN II pro- 
vides only data for the stronger A6584 line. We have 
included the contribution from the A6548 line by as- 
suming a typical ratio A6584/A6548 = 2.94. 

Figures 1, 2 and 3 show the distribution of the 
present dataset compared with the original zones for 
PNe, H II Regions and SNR as proposed by SMB77. 
The present sample requires an expansion of the orig- 
inal axes limits. Figure 1 shows the log [Hq!/[Sii]] 
vs log [Ha/ [Nil]] intensity ratios; Figure 2 the log 
[Ha/[Sii]] vs [Six] AA6717/6731 ratios and Figure 3 
the log [Ha/[Nii]] vs [Sii] AA6717/6731 ratios. In 
Figures 2 and 3 we have added dashed horizontal 
lines that indicate the range for which the [S ii] in- 
tensity ratio is a reliable indicator of the electron 
density (Osterbrok 1974). The corresponding elec- 
tron density scales are included on the right axes of 
these figures. Figure 1 and 4 contain the 613 objects 
in our sample; Figures 2 and 3 and 5 and 6, that 
involve the [Sii] AA6717/6731 ratio contain only 550 
objects for there were 52 objects with only one [S II] 
line listed in the SCGPN II which precluded obtain- 
ing the [S ii] line ratio. Likewise, eleven other objects 
were also discarded in those diagrams for having an 
unlikely large [Six] AA6717/6731 > 2. 

3. DISCUSSION 

It is clear from Figures 1-3 that the limits and 
zones defined for PNe in the original sample used by 
SMB77 fall short of covering substantial data regions 
in these diagrams. The present extended sample 
from the SCGPN II expands the previous location 
of PNe in the diagrams, providing a more accurate 
representation of their zones of influence. 

New statistical limits therefore have now been set 
to represent the location of PNe in the SMB77 di- 
agrams by means of density ellipsoid probabilities. 
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Fig. 1. The database used in this work is plotted in the 
log [Hci/[S ii]] t;s log [Hci/[N ll]] plane. The original zones 
for SNR, H II regions and PNe proposed by SMB77 are 
indicated. 
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Fig. 2. As in Fig. 1, but for the log [Ha/[Sll]] vs [Sll] 
AA6717/6731 plane. The horizontal dashed lines indicate 
the range of values for which the [S ll] AA6717/6731 ratio 
is a valid indicator of the electron density. 

Figures 4, 5 and 6 show the new zones defined by el- 
lipses of probablity 0.85 for the corresponding emis- 
sion line ratios for PNe. The ellipse in the diagrams 
represents both the density contour and confidence 
region of finding an object within this zone with an 
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Fig. 3. As in Fig. 1, but for the log [Ha/[Nll]] vs [Sll] 
AA6717/6731 plane. The horizontal dashed lines indicate 
the range of values for which the [S ll] AA6717/6731 ratio 
is a valid indicator of the electron density. 



Fig. 5. As in Fig. 4, but for the log [Ha/[Sll]] vs [Sll] 
AA6717/6731 plane. The correlation coefficient is -0.58 
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Fig. 4. Diagnostic diagram for the log [Hq/[Sii]] inten- 
sity ratio vs log [Hq/[N ii]] ratio with the new limits de- 
fined by a density ellipse of probability 0.85. The corre- 
lation coefficient is 0.79 

85 % confidence level. The density ellipse is also a 
good graphical indicator of the correlation between 
the two variables involved. The ellipse collapses 
across its minor axis, i.e., becomes thinner, as the 
correlation betwen the two variables increases and 
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Fig. 6. As in Fig. 4, but for the log [Ha/[Nll]] vs [Sll] 
AA6717/6731 plane. The correlation coefficent is -0.37 

approaches unity (±1) and expands or reduces its 
eccentricity as the correlation weakens. The density 
ellipses in Figures 4, 5 & 6 have all been set for a 
0.85 probabilty, however the correlation coefficients 
differ among them, with values of 0.79, -0.58 and 
-0.37, respectively. Thus, the correlation between 
log [Hq!/[Sii]] vs log [Ha/ [Nil]] is strongest whereas 
that between log [Ha/[Nii]] ?;s [Sii] AA6717/6731 is 
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Fig. 7. Histograms of the distribution of the values for 
the corresponding line ratios. Numbers at the top of each 
bar correspond to the percentage value of that column 
in relation to the total sample 



weakest. The correlation is related to the scatter of 
the points in the diagrams and this scatter is a com- 
bination of real dispersion due to different physical 
conditions, as those expected from different evolu- 
tionary stages of a planetary nebula, different condi- 
tions within a nebula (e.g. core, rim or halo) and in 
some cases, also probable observational uncertainties 
or measuring errors from the spectroscopic SCGPN 
II survey. Phillips (2004) has discussed a number 
of factors that potentially affect such samples, such 
as the aperture (slit) size and its location over the 
nebula. 

The samples for each diagram follow the trends 
apparent in the original SMB77 diagrams, though 



now they cover zones where the concentration and 
correlation of the line ratios provide a much closer 
representation of the true location of PNe in these di- 
agnostic diagrams. For example, the ellipses in Fig- 
ures 5 and 6 show a clear tendency towards larger 
electron density values as compared to the corre- 
sponding original areas indicated in Figures 3 and 
4. These values are in good accord with the electron 
densities derived from other, smaller, samples such 
as in StanghcUini and Kalcr (1989). 

The diagrams provide direct and useful informa- 
tion on general physical conditions in PNe. For ex- 
ample, diluted shells with low electron density and 
low to moderate excitation conditions, e.g. 0.5 <log 
[Ha/[Sii]] < 1 (see Fig. 5) will overlap the loca- 
tion of H II regions. These objects are likely large, 
evolved PNe. The likely presence of shock contribu- 
tions in the dominantly photoionized environment of 
PNe is hinted by objects located in Fig. 4 in the log 
[Ha/[Nii]]< and log [Ha/[Sii]] < 0.4 region. Ob- 
jects with high electron density and high Ha emis- 
sivity in Figs. 5 and 6 point towards compact, young 
PNe. For PNe for which limited information is avail- 
able, their location on the diagrams usually provides 
insight into their general (or peculiar) characteristics 
and evolutionary stage. 

Figure 7 shows the histograms of the distribution 
of values for the corresponding line ratios. The num- 
bers at the top of each bar indicate the percentage 
of objects in the sample associated with that value 
in relation to the entire sample. The line ratios of 
PNe are concentrated around a value of 1.2 for log 
[Ha/[Sll]] with a fairly even distribution on either 
side from the peak. For the log [Hq:/[Nii]] line ratio 
its values concentrate around 0.0 but falling rapidly 
to negative values. Finally, the distribution of val- 
ues for the [Sii] AA6717/6731 line ratio is markedly 
asymmetric, peaking sharply at a value of 0.6, with a 
smooth drop towards larger values and an abrupt fall 
on the other side (high density). These values can be 
considered as representative for the most common or 
typical conditions in galactic PNe. 

4. CONCLUSIONS 

Values of Ha/ [Nil], [Sii] AA6717/6731 and 
Ha/[Sii] for 613 planetary nebulae have been de- 
rived from the Strasbourg Catalogue of Planetary 
Nebulae, part II (Acker et al. 1992) and plotted in 
the Sabbadin, Minello & Bianchini (1977) diagnos- 
tics diagrams. The current sample substantially ex- 
pands the sample used in the original diagrams and 
allows a redefinition of the parameter space of PNe 
in the log [Ha/[S ii]] vs log [Ha/[N ii]], log [Ha/[S ii]] 
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vs [Sii] AA6717/6731 and log [Ha/[Nii]] vs [Sii] 
AA6717/6731 planes. Accordingly, new statistical 
limits have been set for the location of PNe in these 
diagrams using density ellipses of probability 0.85. 
We find a good, 0.79 correlation for the log [Hq;/[S ll]] 
vs log [Ha/[Nll]] intensity ratio; a modest -0.58 cor- 
relation for the log [Ha/ [Sii]] vs log [Ha/ [Nil]] in- 
tensity ratio and a poor, -0.37 correlation in the case 
of log [Ha/[Nii]] vs [Sii] AA6717/6731. Prom our 
study wc obtain the statistical distribution of emis- 
sion line ratios expected in PNe. These upgraded 
diagnostic diagrams provide a better representation 
of the ranges of physical conditions indicated from 
the emission line ratios and allows the identification 
and study of groups or individual PNe of particular 
interest in specific regions of the diagrams. 

Finally, in relation to the parameter space for the 
H II and SNR regions, the data sources in SMB77 for 
these two cases were restricted to a limited number 
of objects, as in the original case for PNe. For the 
HII regions they used the compilation made by Kaler 
(1976). For the supernovae the data is based on 17 
galactic remnants. It would be clearly useful to re- 
vise and update the parameter space for the H II and 
SNR regions to assess the reliability of their distribu- 
tions in the revised diagrams. This is presently out 
of the scope of the present paper but worth pursuing 
as a natural extension of this work. 

Templates in poscript or pdf format of the revised 
diagrams. Figures 4, 5 and 6, are available from J. 
A. Lopez. 
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